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Abstract 

The bound eigenfunctions and spectrum of a Dirac hydrogen atom are found taking 
advantage of the SU(1, 1) Lie àlgebra in which the radial part of the problem can be 
expressed. For defining the àlgebra we need to add to the description an additional 
angular variable playing essentially the role of a phase. The operators spanning the 
àlgebra are used for defining ladder operators for the radial eigenfunctions of the 
relativistic hydrogen atom and for evaluating its energy spectrum. The status of the 
Johnson-Lippman operator in this àlgebra is also investigated. 

Key words: relativistic hydrogen atom, ladder operators, SU(1, 1) Lie àlgebra 
PACS: 33.10.C, ll.10.Qr 



The bound solutions of the hydrogen atom are of great importance in both 
classical and quantum mechanics and so is the search for new ways of solving or 
using such problem [1,2,3,4,5,6,7,8,9,10]. The purpose of this paper is to discuss 
an algebraic solution for the bounded eigenstates of the relativistic hydrogen 
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atom based on the SU (1,1) properties of the extended radial part of the 
problem [11,13]. The àlgebra is spanned by two operators Sj, i = 1, 2 in terms 
of which such part of the problem can be expressed, plus an additional — hence, 
extended — phase shifting operator, S 3 , whose only effect is to change the 
relative phases of the eigenfunctions of the àlgebra participating in its radial 
eigenfunctions; though, for maintaining the eigenfunctions of the hydrogen 
atom, it is found necessary to project our results to vanishing additional phase 
after applying the operator S 3 . That notwithstanding, the set of operators 
allow a very direct solution of the relativistic hydrogen problem. Moreover, 
the solution exhibits that the extended radial system possesses a SU (1,1) 
symmetry 1 — the original problem possessing an 5*0(2, 1) one [14]. We further 
exhibit how the eigenfunctions of the àlgebra, i. e. the radial eigenfunctions 
of the relativistic hydrogen atom, can be related to the Laguerre polynomials 
of noninteger indices and to the Sonine polynomials [15]. These results can be 
of importance in the study of squeezed and coherent states of the relativistic 
hydrogen atom, as well as in diverse quantum òptics applications [6,16,17,18]. 

That the relativistic Coulomb problem could be solved exactly and that the 
energies so calculated showed an excel·lent agreement with the experiments 
was indeed one of the reasons for the ràpid acceptance of Dirac relativistic 
quantum mechanics. Almost from the start operational methods were applied 
to the hydrogen atom, as in [19,20,21]. In due course other methods were 
developed, for example in [22,23], and in Biedenharn's solution [24] in which 
advantage was taken of the constant operator B introduced by Johnson and 
Lippmann [24,25,26] — shown in equation (42) below. 

The Hamiltonian of the relativistic hydrogen atom is 

H D = ca-p + l3c 2 --, (1) 

r 

using atòmic units (H — m — e — 1), c is the speed of light, a and f3 are the 
Dirac matrices in the Dirac representat ion [27,28], 




where the l's and O's stand, respectively, for 2 x 2 unit and zero matrices and 
the a is the operator composed by the three Pauli matrices a = (<7 x ,a y ,a z ). 
The method to be discussed in this paper is different from all the aforemen- 
tioned, as it deals exclusively with the radial part of the problem. Besides, to 
obtain a system of equations invariant under SU (1, 1), we need to extend the 

1 symmetry that can, by suitably changing the definition of the Sj — hence the 
conmutation relations — be recasted as a SU(2) one[7,ll,12] 
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configuration space by adding an extra angular variable ip (Cf. equation (17) 
below). 

Since the Hamiltonian (1) is invariant under rotations, we look for simultane- 
ous eigenfunctions of Hp, |J| 2 and J z , where J = L + s is the total angular 
momentum, and s = |S is the spin. Parity is a good quantum number in our 
problem, but instead of parity, we use the parity related quantum number e, 
defined as 

'i in = j + |, 

(3) 



we can easily also show that the eigenvalues of the operator K = (3(1 + E • L), 
are given by —e(j + 1/2), so, in a way, e is a more useful quantity than parity. 

The lower components of the wave function should have an opposite parity 
than the upper ones, and as parity goes as (—1)', where / is the orbital an- 
gular momentum quantum number, if we define V — j — e/2, we can cast the 
eigensolutions of (1) as 

1 / F(r)yT=i-e,2Í0A) \ 
*(r,M) = - ■ (4) 

r \iG(r)yf =v+e/2 (eA)) 



where y™ = i± e / 2 (@, <P) are spinorial harmonies [29]. The form proposed for the 
eigenfunctions solves the angular part of the hydrogen atom. 

Let us define the quantity k = y/c 4 — E 2 , which is positive definite since we 
are interested in bound states. Furthermore, let us also define 

/ 1 \ I c 2 — E 

( = a F =l/c, Tj = e^j + -j, v = y ^q^, (5) 



where cíf — 1/137 is the fine structure constant. With the proposed solution 
(4), we have to deal only with the radial equation for the relativistic hydrogen 
atom [30] 



H 



Dr 



a r 



Pr ~ Í— (3 

r 



+ c 2 (3 



(6) 



where 

a r = -a-r, p r = —- I 1 + r- 1 I . (7) 
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Using the radial part of the solution (4), definitions (5) and (7), equation (6), 
and using p = kr, we can obtain the pair of equations [8] 

and 

The above equations, (8) and (9), can be rewritten using a set of three op- 
erators whose commutation relations define a SU(1, 1) àlgebra, as we exhibit 
in what follows. Introducing the variable x through the relation p = e x , so 
x G (—00,00), and the functions ip±(x) through 



F(p(x)) = Vc 2 + E [ÍJ_(x) + i) + {x)\ , 
G{p{x)) = Vc 2 -E - ^+(x)] , 



(10) 
(11) 



we obtain 
d 

dx 



CE 



Vc 4 - E 2 



tp + (x) 



C 



;i2) 



and 



eix 



Vc 4 - £ 2 



c 



Vc 4 - £ 2 



+ r i U+(x). (13) 



This first-order system can be uncoupled to yield 

Mx) = (r^-C 2 -\)Mx), 



& 2+ ^-^~\ 



(14) 



and 



dx 2 4 



V>-(x) 



(jj-ç 2 -\)l>4*), (is) 



where we have defined 



Ve 4 - -E 2 



+ 1/2. 



(16) 
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The seemingly odd —1/4 term in (14) and (15), has to be included in order 
to close the àlgebra that solves the problem. We pinpoint that the change 
of variable from p to x is, strictly speaking, not necessary for any of the 
calculations that follow; the change, however, simplifies the appearance of 
some of the equations. Notice also that we can regard the equations (14) and 
(15) as eigenvalue equations with a common and predetermined eigenvalue 
£ = — ( 2 — j = j(j + 1) — ( 2 . Since the minimum value of j is 1/2, then 
£ > 0. In cases where the atòmic number, Z, is different from 1 — the value 
we are assuming in this article — then £ > for Z — 2, • • • up to 118. If Z 
exceeds this value then £ could become imaginary and the solutions would 
then exhibit irregular behaviour [26]. 

To construct the Lie àlgebra for our system, let us introduce an extra variable, 
(f G [0,27r], through the operator 

(17 > 

and then define 



— ± 



— e x i— + - 

dx dip 2 



^k^T^ + í > ( 18 ) 



which depend both on (p and on the transformed 'radial' variable x. These 
operators satisfy the àlgebra 

[S 3 ,S ± ]=±S ± , [~ + ,~_] = -2~ 3 . (19) 



Let us now introduce the operators Si and S 2 

2i = Í(S + + S_), S 2 = 1(S + -S_), (20) 



which together with S 3 satisfy the SU(1, 1) Lie àlgebra [4,13], 

[Si, S 2 ] = -iE 3 , [S 2 , S 3 ] = iSi, [S 3 , S^ = iE 2 ; (21) 

whose Casimir operator is [13] 

S c = -El -E 2 2 + E 2 3 = d 2 /dx 2 - e 2x - 2ie x d/d<p - 1/4. (22) 



Note that we do not have a linear term in d/dx in the Casimir operator 
precisely because we included the 1/2-term in the definition of the S ± . It is 
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also worth noting that the realization of the S77(l, 1) àlgebra given in the 
previous equations is in terms of two variables (another example is described 
in [31]) and not in terms of a single one as it is more usual [32,33]. 

To find the simultaneous eigenfunctions and eigenvalues of the operators S c 
and S 3 , let us change to the variable x, using p = e x , and write the eigen- 
functions as Vç(x,(p); where /i and £ stand for the eigenvalues of S 3 and S c , 
respectively, 

Z 3 V£(x,ip) = »V£(x,<p) 

Z c V£(x,<p) = tV£(x,<p). (23) 

Using equation (17), we find 

V£(x,(p) = e i '«'Pg(x). (24) 

It is clear that f plays the role of a phase, and 

Mx) = PÇ(x), j,-(x) = PT\x). (25) 

From equation (19), we can also conclude that the S± are ladder operators 
for the problem, they change /i to /i ± 1, so they are easily shown to satisfy 

S ± V^x,ip) <xV^\x,ip). (26) 

To establish the properties of the S-operators we need to define the inner 
product 

2-7T OO 

<0,^>= í í (f>*(<p,x)il>(<p,x)dx, (27) 



2tt 

-oc 



where we are assuming that <f>(ip,x) and ip(ip,x) are periòdic functions over 
the interval ip G [0, 27r] and that they vanish as x — > ±oo — so as to make them 
square-integrable. Using (27), we can prové that S 1; S 2 , and S 3 are hermitian, 
"i = Cj, i = 1,2, 3. 

To see that the problem of finding the radial eigenfunctions of the extended 
radial part of the relativistic hydrogen atom is equivalent to finding a repre- 
sentation of the SU(1, 1) àlgebra, we introduce a complete orthogonal basis 
of simultaneous eigenfunctions for S c and S 3 , namely V^(x,ip) = |£/i >, we 
require that < A* >= ^h,h'j using the inner product (27). 
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The operator E 2 = S 2 + S 2 , + H 2 is obviously definite positive; furthermore, as 
E 2 = 2~ 2 - S c , (28) 

we must have 2/i 2 > £, this in turn means that is bounded by below, as 
a consequence there is a minimum value for it, say |/i| m j n . We have then two 
choices, either we choose /i > or we choose fi < 0. In the first case /i itself is 
bounded by below; in the second, ji is bounded by above. 

Let us take first /i > and define À = /i m j„. Let us call À = |/i| m i„ then 
2_|£À >= 0, or equivalently S+S_|£A >= 0. Since S + S_ = — S c + S| — S 3 , 
this means 

_£ + A 2 -À = or £ = A(A-1). (29) 
As À has to be positive, we directly obtain that H m { n = A = s + |, where we 



have defined s = +yrf — ( 2 = +y (j + 1/2) — a 2 ?, so s is real. 

Some features of this algebraic approach are worth noting. Fist, notice that the 
eigenvalue £ = A(A — 1) is not determined by the àlgebra itself, but just by the 
angular symmetry of the problem. Second, that this eigenvalue is not necessar- 
ily an integer or half-integer number as it should be in the related SU (2) case. 
Nevertheless, it is interesting to point out that in the non-relativistic limit À 
becomes approximately half-integer since in such a limit we can neglect the 
fine structure constant «p, so s ~ (j + 1/2). In this case the eigenvalues of H c 
become 1(1+ 1) but the structure of the àlgebra 577(1, 1) remains unchanged. 
The remarkable point here is that the representat ion of the àlgebra and the 
energy spectrum become precisely those of the non-relativistic hydrogen atom 
[45]/ 

With the value of À fixed, we proceed to obtain the representations of £77(1, 1). 
We found it convenient to relabel the kets >, using £ = A(A — 1), as 
\Xfj, >. So, starting with |AA >, we obtain an infinite — since \i is bounded 
only by below — series of states |A A + 1 >, |A A + 2 >,•••, by aplying 1, 2, • • • 
times H + to the ground state |AA >. Therefore, /i itself takes the vàlues 
[i = A, A + 1, A + 2, • • • = A + n, with n = 1, 2, • • • 

To find the energy spectrum, we just need to solve, from definition (16), E in 
terms of /x; thus the energy spectrum of the hydrogen atom is 



E = c 



1+ " 



(/i - 1/2)= 



-1/2 

= c 2 



1 + 



(/i - 1/2)= 



-1/2 

(30) 



where /j, = X + n, or /i — 1/2 = s + n, s — \J (j + 1/2) 2 — a 2 F , n — 0, 1, 2, 
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Let us write the above reult in a more familiar form. We define the principal 
quantum number N — 1, 2, 3 • • • and the auxiliary quantity ej as follows 

N = j + l/2 + n, 

€j=N-s-n = j + 1/2- s. (31) 

From this definition we have that p — 1/2 = s + n = N — €j. When we 
introduce this last result in the energy equation (30), we obtain precisely the 
Dirac energy spectrum for the Hydrogen atom,[21,27,28, 30,35]. 

To obtain the eigenfunctions from the representation, let us write 

Z±\\p >= C±\\p± 1 > . (32) 

We can evaluate these constants from < À //|S + S_|A p >= — À(A — 1) + p(p — 
1) = C^_ X C~ and from (C^)* = C~\ if we further assume these constants 

to be real, we get = + ^J p(p ± 1) — À(À — 1). 

The ground state wave function is obtained from S_|A A >= 0, that is 

p£( x ) = N x e sx exp{-e x ) = N x e^' 1 ' 2 ^ exp(-e x ) = N x p s e~ p , (33) 



where N\ = 2 / A 1/2> is the normalization constant calculated using the in- 

v /r(2A-i) to 

ner product (27) and T(2À — 1) stands for the Euler-gamma function. From 
equation (25), we find that for the lowest eigenvalue, ip+(x) = P x (x), and 
if)-(x) = 0. In terms of the radial variable p, the ground state solutions are 
[equation (10)] 



F(p) = N * 



^ 2 + E) .p s e^, G{p) = -N, 



\2c\\-l/2) F ' yy '~ ^2c 2 (A-l/2) 



(C2 - E) p s e~ p . (34) 



The normalization constant for this state, N = 1/y 2c 2 (À — 1/2), different 
from N\, is obtained from the condition 

oo 

/ (\F(p)\ 2 + \G(p)\ 2 )dp=l. (35) 



The rest of the eigenfunctions for the positive set of eigenvalues p are obtained 
by applying successively S + to the ground state |A A >. Notice that the value 
of À is determined by the total angular momentum and is different for each 
value of j. 
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The representation presented here is called the two-mode bosonic representa- 
tion, because the operators (18) can be constructed using two bosonic creation 
operators [46] . This can be easily seen if we define the two mode bosonic oper- 
ator operator A = l/2(à^à — b'b+1), where and ¥ are the creation operators 
of (some) bosonic modes and à and b are the anhilation ones. If we put \m,n> 
as the simultaneous eigenstates of the operators à^à and ¥b, respectively, we 
can see that our representation is given by \X p >= \n + 2À — 1 ,n >, with 
p — X + n, n — 0, 1, 2 • • •. In the non relativistic case À becomes an integer (the 
principal quantum number becomes an integer), because À = l + l, [45] and the 
representation becomes what Bargmann calls minimal M, (X has a minimal 
value), where À takes the vàlues 1/2, 1, 3/2, • • -[32]. In the non- relativistic 
case the semi-integer representat ions are excluded [45]. It is interesting that in 
the relativistic case we use a generalization to a non-integer, nor semi-integer 
index, of the minimal representation. 

What happens when we consider the negative set of eigenvalues of S 3 (/i < 0)? 
In that case |//| m j n = A is the largest eigenvalue, instead of the smallest one. 
Consequently, we obtain the eigenfunctions by successive applications of S_ 
to the ket |AA >. Starting with this state and repeating the procedure we 
follow in the p> case we find that 

|A A >~ e- iX ^e sx e eX = e^p'e? (36) 

Which diverges as p — > oo. This behavior makes the negative eigenvalues So- 
lutions not square integrable, and we must have to discard them if we want 
to describe physically realizable states, keeping only the positive set of eigen- 
values of jl. . 

As a last point, we note that in [7] it is shown that all the solutions are 
polynomials with weight factor W(p) = p x ~ 1 l 2 e~ p '. This weight factor assures 
that the behavior of the big and the small components of the spinor are regular 
both at the origin as well as p — > oo. 

We have thus constructed an SU (1,1) àlgebra for the relativistic hydrogen 
atom by introducing the Hermitian operators Sj. The eigenstates of the prob- 
lem are labelled by numbers p, which are neither integers nor half-integers. 
Let us remark that we were forced to introduce the new variable ip in order to 
close the àlgebra. In terms of the solutions of the Dirac equation, ip just plays 
the role of an extra phase. The functions, í/j±, defming the eigensolutions of the 
hydrogen atom [Eqs. (10) and (11)] are, in strict, obtained by the projection 
to zero additional phase of the eigenfunctions of the àlgebra, namely 

Mx) = PÇ(x) = V£(x,y% =0 , ^_(x) = P<TV) = Vf-W)U=o- (37) 
What happens to the radial solutions F(p) and G(p) when we apply the op- 
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erators S c and S 3 to the eigensolutions V^(x,ip)l Given that both V^(x,ip) 
and Vç t ~ 1 (x, ip) are eigenfunctions of S c , the energy spectrum and eigenfunc- 
tions remain essentially unchanged, excepting for an irrelevant global constant 
À(À — 1). But, when we perform a S 3 rotation things are different; in this case 
we change F(x) to 

F(x) [e^-^-ix) + e^ip+ix)} = é^\e~^i)_{x) + ip + (x)}, (38) 

since F(x) ~ [ip^(x) + ip + (x)]. However, the energy spectrum does not change 
since both V^(x, ip) and S^V^x, ip) have the same eigenvalue [i. Nevertheless, 
as we can see from equation (38), though the phase e* w does not play any 
observable role, the term e~ tlp changes the relative phase between the ip+(x) 
and the ip_(x) and so it changes the radial function F(p) itself. This means 
that the rotated eigenfunction would no longer be an eigenstate of the radial 
hydrogen atom Hamiltonian unless we first project the solution to vanishing 
phase: <p — > 0. We could say that there should be a sort of superselection rule 
between the functions ip-(x) and ip + (x), this interpretation, apart from pro- 
viding a way out of that minor problem, would also confirm that the operator 
S3 cannot be an observable[36,37]. All these peculiar features of our solution 
method might have some bearings for the properties of certain other systems 
[38]. 

We remark that the eigenfunctions ip+{x) = P\(x) and i^-{x) = P^ ^^x) are 
related to the Laguerre polynomials of non integer index [21], Lf l s (2e x ) and 
L^e*), as [30] 

P»(x) = ae sx e~ eX Ll s (2e x ), PT\*) = be sx Ll^e*) (39) 
where fj, — n + s + 1/2, and the constants a and b should satisfy 

a(Tj + s- (W 1 + n) + b(n + 2s) = 0, 

b(Tj -s + -n)-an = 0. (40) 

In fact, if we solve for a in the second equation and substitute it back in 
the first of (40), using also that s + n = ji — 1/2, we again obtain precisely 
the energy spectrum of equation (30). Notice that these non-integer Laguerre 
polynomials are related to the hypergeometric function 1 F 1 (— n, a + 1; x) and 
to the Sonine polynomials T^ n \p) through the relations [15] 

Ta . . T(a + n + 1) _, . 1 . 

L ^= nir(a + l) W -»>° + 1 ·>P) 

= (-l) B r(a + n + l)TW(p). (41) 
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Let us close the paper with some considerations about the Johnson and Lipp- 
mann operator, B, defined as 



B = ^-K l5 (H D -f3c 2 )--Z-r, 



(42) 



where 75 = ia1a.2a.3-, and K = /3(1 + E-L). This operator was introduced many 
years ago to study the properties of the Dirac Coulomb Hamiltonian [25,24]. A 
direct calculation shows that the Hamiltonian of the relativistic hydrogen atom 
commutes with the Johnson and Lippmann operator, which can be regarded 
as the relativistic counterpart of the Laplace-Runge-Lenz vector of classical 
mechanics [34]. 

To calculate the action of B over the state functions (4), we need to use the 
identities (a ■ r)y^ l±e / 2 (d, 4>) — — 3^=^/2 4Ò'i a fter some àlgebra, we found 
that 



B*(r,0,0) 



dr 



-y}m 



(M) \ 



\ -ít.í í -dG i TjG ÇF\ , CiG l -y 
\ [ r V dr ~ r r ) ' r J ^ ; 



y j=l+e/2 



(43) 



Here we see that the efect of this operator over the angular part interchanges 
the coupling of l with V but the total angular momentum quantum number j 
remains unchanged. This point explains one of the degenerades of the wave 
functions in the relativistic Coulomb problem, since the energy of a particular 
state depends on iV and j and not on the coupling of the orbital angular 
momentum with spin. 

We calculate now the effect of B on the wave function using the equations 
of motion (8) and (9), changing to the variable p = kr and writing only the 
effect over the radial part of the function 



B 



I ( F (P) \ 
p\iG(p)J 



c 2 p 



ÍTj(c 2 -E) 



( F{P) 
\iG{p) 



(44) 



As expected, the Johnson and Lippmann operator on shell is no longer a 
differential operator. Unfortunately, the operator is not diagonal in terms of 
F and G; nevertheless, using the functions ip±, we can write the above result 
as 



B 



I ( F (P) \ 
p\iG(p)J 



1 V^Te [(c - f )V + + (C + f )V> 
p { V^e \-(c - ^)^ + + (C + *£)1>- 



(45) 



again we found a non diagonal result, nevertheless, from equation (25), we 
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know that at least the ground state |AA > correspond to a diagonal case, 
because for that state = 0. Are there any other state in which this feature 
occurs? In order to obtain a diagonal state, we need that one of the two choices 



which is never null excepting again for the ground state |ÀÀ >. But this is 
the ground state of the reprès entation of the SU(1, 1) àlgebra and not of the 
hydrogen atom. In terms of the wave functions of the relativistic Coulomb 
Hamiltonian, the ground state of the àlgebra corresponds to the infinite set of 
states. Namely, (N = 1 and j = 0), (N = 2 and j = 1), (N = 3 and j = 2) and 
so on. In other words, using spectroscopy notation iVj, we find the interesting 
result that the Johnson and Lippmann operator takes a diagonal form when 
written in terms of ip±, only for the states ls, 2p, 3d . . . 

In conclusion, we have solved exactly the Dirac relativistic hydrogen atom 
using an algebraic technique based on the SU (1,1) Lie àlgebra in terms of 
which the radial part of the problem can be expressed. With its help we 
calculated the explicit form of its ground state wave function. We have also 
derived radial ladder operators for the system and pointed out the way for 
obtaining all of its radial eigenfunctions using them. Let us also pinpoint that 
in the limit ( — > the àlgebra (19) [or (21)] trivializes, /i — > 1, and the 
spectrum becomes continuous, as corresponds to a free Dirac particle. 

The method discussed in this paper, apart from the algebraic solution of the 
relativistic hydrogen atom, could be of interest in itself as it requires extending 
the system before the solution and then projecting the results back to zero 
phase tp = — although this last step is not always necessary. It could be an 
interesting problem to account for all relativistic problems that may be solved 
by a similar extension and to find their corresponding extended symmetry 
algebras, as well as to investigate their generalization to D > 3 spatial dimen- 
sions [39]. The algebraic method of solution and the expression of the bound 
state eigenfunctions in terms of the ladder operators S± presented here may 
offer some simplifications in relativistic atòmic physics calculations [40], can 
be useful in diverse applications in quantum òptics [6,17] and it may, possibly, 
be applied to the study of ionisation states [41] by using non-discrete repre- 
sentations of 577(1, 1); at least below the pair-production treshold, or in the 
non-relativistic limit [42,43,44], where a description of the hydrogen atom in 
terms of such an SU(1, 1) àlgebra is also possible [45]. 

As we already said before, other application of the method described here 





(46) 
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concerns the non-relativistic problem [42,43,44], where a description of the 
hydrogen atom in terms of such an àlgebra is also possible. The point here is 
that the 577(1, 1) àlgebra remains exactly the same also in the non-relativistic 
limit, and only the relationship of /i with the energy, the value of s itself, and 
the fact that we have to deal with only one radial function instead of two, are 
the only important points we need to change [45]. 
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